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ABSTRACT 

General  expressions  are  developed  for  the  fluctiiation  in  density 
of  electrons,  ions,  and  charge  in  a  plasma  in  thermal  equilibrium  in 
an  extemed  magnetic  field  taking  only  Coulomb  interaction  into  account. 
The  spectral  distribution  of  the  spatial  Fourier  components  of  these 
fluctuations  is  derived  from  basic  principles. 

The  fluctuations  in  electron  density  are  discussed  in  some  detail, 
and  spectra  are  confuted  under  conditions  which  are  thought  to  prevail 
in  the  outer  ionosphere.  Frequency  spectra  of  general  validity  are  com¬ 
puted  for  electron  density  fluctuations  along  the  magnetic  field.  It 
is  shown  by  means  of  examples  that  the  frequency  spectra  under  ionospheric 
conditions  sure  not  much  Influenced  by  the  magnetic  field  except  for  den¬ 
sity  fluctuations  fairly  close  to  perpendicularity  to  the  magnetic  field. 

Applications  to  Incoherent  backscatterlng  sure  discussed,  smd  it  is 
shown  that,  under  suitable  conditions,  backscatter  techniques  can  give 
valuable  information  about  electron  density,  temperatvure  and  consti¬ 
tuents  of  the  ionosphere. 
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I.  mTRODDCTiar 


Much  of  the  current  Internet  in  the  fXueti)knti<m  of  electron  den¬ 
sity  in  e  plasm  In  themodynanlc  equilibrium  stems  from  an  original 
suggestion  by  V.  E.  Gordon  (19^)  that  such  fluctuations  might  be  de¬ 
tectable  in  the  upper  ionosphere  in  backscatter  observations  with 
powerful  radar  systems.  Gordon  predicted  the  toteG.  scattering  cross 
section  by  adding  up  the  cross  sections  of  the  individual  electrons, 
and  he  estimated  the  spectraG.  broadening  by  ascribing  to  the  scatter¬ 
ing  elements  the  thermal  velocity  of  individual  electrons. 

Hhen  Bowles  observed  this  type  of  scattering  19^9)  li 

appeared  that  the  total  power  was  of  the  ri£^t  order  of  magnitude,  but 
that  the  spectrum  was  noidiere  near  as  wide  as  expected.  This  has  been 
confirmed  by  recent  observations  by  Plneo,  et  al,  (i960) .  Subsequently 
sevooral  authors  have  desoribed  theories  to  explain  this  dlscrepney 
and  hanre  eoae  ^  vllh  essentially  similar  results  (Tejer  i960,  Salpeter 
1960a,  1960b,  Oou^^rty  and  I^ley  1960a,  1960b) .  The  reason  for  the 
narrow  spectrum  is  the  interaction  between  the  electrons  and  the  lone. 
This  interaction  causes  the  spectral  width  to  be  that  which  would  have 
resulted  if  the  electrons  had  had  the  mass  of  the  ions  and  had  been 
completely  free  to  move,  very  roughly  stated.  The  theories  referred 
to  show  that  the  actual  spectral  distribution  is  considerably  more 
complicated  if  all  the  details  are  taken  into  account.  Only  Douc^rty 
and  Farley  have  taken  the  effect  of  the  geoomignetic  field  properly 
into  account. 

The  present  work  is  concerned  with  the  calculation  of  the  electron 
density  fluctuations  in  the  gmeral  case  of  ions  with  an  arbitrary 
number  of  positive  chsurges  and  in  the  presence  of  an  external  magaetic 
field.  The  effect  of  two-body  collisions  is  taken  approxlmtely  into 
account  by  introducing  a  relaxation  term  in  the  Boltsaann  transfer 
equation.  The  derivation  does  not  depend  on  Byqulst's  theorem  as  does 
that  of  Douid^rty  and  Farley,  but  their  result  is  obtained  in  the  case 
of  singly  charged  low  and  no  collisions.  When  the  magnetic  field  is 
neglected,  but  the  multiple  ion  charge  retained,  the  result  of  Salpeter 
(1960b)  is  recovered. 
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i  The  calculation  of  the  total  fluctuation  Is  first  carried  out  along 

lines  corresponding  to  those  used  by  Pines  and  Bobm  (199^) .  These 
•  total  fluctuations  are  shown  to  be  lndm>end0nt  of  an  external  nagnetlc 

field  and  Independent  of  t«D»body  eolllelotti«  Tins  Variations  of 
the  density  fluctuations— -or  the  spectrum— are  then  discussed  by 
means  of  the  Boltzmann  equation.  The  solution  of  this  equation  fol¬ 
lows  very  closely  that  of  Bernstein  (1958)  who  investigated  waves 
in  a  plasma  In  a  magnetic  field.  We  essentially  have  only  to  supply 
the  appropriate  Initial  conditions  to  the  solutions  given  by  Bernstein. 
It  turns  out  that  the  fluctuations  In  ion  density  and  chetrge  density 
and  their  spectral  distribution  may  be  obtained  with  very  little  extra 
work.  These  results  are  therefore  given  In  addition  to  the  electron 
density  spectra.  Tliey  are  of  no  direct  interest  to  the  radar  back- 
scatter  observations j  but  may  be  of  use  in  other  problems,  such  as 
determining  friction  and  diffusion  coefficients  of  the  Fokker-Planck 
eqxiatlon  for  the  plasma  (Hubbard  i960). 

Numerical  calculations  of  the  spectral  distribution  of  the  elec- 
tron  density  fluctuations  are  made  for  a  set  of  paraaieterB  covering  a 
large  range  of  temperatiures,  plasma  frequencies  and  radar  frequencies 
appropriate  to  ionospheric  conditions.  The  results  are  given  as  sets 
of  curves  which  can  be  conveniently  used  for  interpreting  backscatter 
observations . 

The  extension  of  the  theory  to  a  larger  nualber  of  Ionic  consti¬ 
tuents  Is  straightforward  in  principle,  but  is  not  carried  out  here. 

II.  RELATIQIB  BEEKEBN  PLABNA  FLUCTIMTIQNB  AID  SCATTBlIHa  FR0PBRTIB8 

Before  we  steurt  discussing  plasma  properties  we  must  relate  the 
scattered  energy  to  the  electron  density  fluctuations.  Por  weak 
scattering,  l.e.,  when  the  Bom  approximation  can  be  used,  one  ob¬ 
tains  for  the  scattering  cross  section  per  unit  solid  angle,  per  unit 
incident  power  density  and  per  unit  scattering  volume 

a  =  OgV  <|n(l)|^>^^  (1) 
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\^ere  Is  the  scattering  cross  section  per  unit  solid  angle  and  per 
unit  Incident  power  density  of  a  single  electron^  and  where  ^  Is  the 
difference  between  the  wave  vectors  of  the  lncld«it  and  the  scattered 
wave.  Note  that  c  depends  both  on  the  polarization  and  the  wave  vec- 
tor  difference  K.  The  q\aantlty  n(k)  is  defined  as 

n(^)  »  ^  5  exp(lto) 

V 

where  n(r)  is  the  number  density  of  electrons.  If,  In  addition,  one 
requires  the  power  spectrum  of  the  scattered  energy,  one  must  study 
the  scattering  cross  section 

=  (^e  V<|n5c,cD)l\v  ^2) 

where  o)  -Is  the  difference  of  the  euigular  frequency  of  the  Incident 

wave  and  the  scattered  power  under  consideration.  The  quantity 

<|n(5,<o)|^  Is  singly  the  power  spectrum  of  the  spatial  Fourier  cooH' 
av 

I»nent  of  wave  vector  k.  The  two  scattering  cross  sections  sure  re¬ 
lated  throu^  the  equation 

+00 

a  =  j  do,  (3) 

-00 

Having  now  established  the  relationships  between  the  scattering  pro¬ 
perties  of  the  plasma  and  the  density  fluctuations  we  can  turn  to  the 
study  of  the  latter. 

III.  TOTAL  FLUCTUATIONS  IN  THE  PLASMA 

Assume  that  the  plasma  Is  neutral  and  that  the  average  number 

density  of  positive  Ions  Is  and  the  average  number  density  of 

electrons  Is  n  .  The  number  of  charges  on  the  Ions  Is  Z  >  n_/N_. 

o  o  o 

The  densities  may  be  eiqpressed  as 


V 

and 

N(r) 

-I 

1«1 

n  V 

0 

IT' 

6(r  -  tj) 

n(r) 

-I 

i«l 

Hr  - 

(h) 


Here  V  is  the  volvnne  of  a  large  periodicity  cube  with  sides  L,  and 
and  r^^  are  the  positions  of  the  ions  and  the  electrons.  The 
charge  density  then  becomes 


p(r)  =  e  tZ»(r)  -  n{r)] 


(5) 


and  the  spatial  Fourier  conrponent  becomes 

P®  -  e  tZIl(k)  -  n(^)]  (6) 

with  N(k)  =  ^  ^  d(r)  N(r)  exp  (iSr) 

V 

and  with  a  similar  expression  for  n(^).  The  wave  vector  k  is  given 
by  k  ■  where  the  i's  are  integers  between  -«  and 

•»«. 

The  intersiction  between  the  various  cheurged  particles  may  now  be 
taken  into  account  through  an  electric  field  ^(r)  which  may  also  be 
expanded  within  the  periodicity  cube  V  *  in  a  Fourier  series. 

By  far  the  most  in^rtant  interaction  in  a  nonrelativistlc  plasma  is 
through  the  Coulomb  forces.  This  means  that  we  let  the  velocity  of 
interaction  become  infinite  so  that  the  electric  field  may  be  derived 
frcxn  a  scalar  potential  idiich  can  be  found  from  Poisson's  equation. 

We  therefore  obtain 

- - 5  P  (^) 

^o^ 


1*  - 


(8) 


t 


f 


vlth  M  dielectric  constant  In  vacuo.  This  vlU  be  a  good  apiproxl- 
mtlon  as  long  as  the  thermal  energy  of  the  electrons  Is  considerably 
swller  than  the  relatlTlstlc  rest  energy  of  the  electrons,  l.e., 

KT 

—5  «  1. 

me 

The  total  energy  of  the  plasma,  inclading  the  self  energy  of  the  indi¬ 
vidual  particles  and  their  kinetic  energies,  then  becomes 

N  V  n  V 

o  o 

i-1  i=l  V 

(9) 

where  M  is  the  mass  of  an  ion  and  m  that  of  an  electron.  Using 
Parceval's  theorem,  the  last  part  of  this  may  be  eigpressed  in  terms 
of  Fourier  conponents  as  follows: 


We  note  in  particular  that  this  expression  will  remain  the  same  irtiether 
there  is  an  external  magnetic  field  present  or  not,  and  it  will  not  be 
altered  by  the  presence  of  neutral  particles  with  which  the  ions  and 
the  electrons  may  collide. 

If  K(^  and  n(r)  were  continuous  functions  of  position,  we 


know  from  Information  theory  that  in  order  to  determine  the  Fourier 
coaip(»ents  vqp  to  t  >  2x(i^,i2,i^)/L  the  number  of  saapling  points 
req]alred  in  space  is  very  close  to  (Brlllouln  .1956) .  If, 

therefore,  means  that  many  particles  must 

contribute  to  each  sampled  value.  Let  these  sampled  values  be  denoted 


by  ... 


for  ions  and  n^ 


*'8(i  It)  electrons. 


-  5  - 


It  it  of  Itttereat  to  knov  bov  these  saapled  values,  or  occupation 
nuBA>ers,  are  related  to  the  actual'discontlnxioue  functions  lf(f)  and 
n(r).  Consider  wave  nuBdbers  ^  where  In.  I,  In. I,  and  |  nJ 

are  saaller  than  or  equal  to  am  iy  respectively.  The  nusflMr 

of  saaqpling  points  required  along  the  three  axes  is  then  2i^  1, 

dig  1  Again,  from  information  theory,  it  follows  that 

the  sampled  values  (occupation  numbers)  may  be  obtained  from  N(r)  and 
n(r)  by  integration  over  the  periodicity  cube  with  the  following 
weighting  factor: 


f(r-r  )  ■  n 

V"2>“3  iil 


0  4  n 

Y.( 

z'  m.L  N 

0XXA 

^  1  2ij^+l  J 

(fiij^+1)  sin 

Vx 

(11) 


If  we  denote  the  sampled  values  at  r 


K  ”2  ^  I 


n(r_  _  «  ),  it  can  be  shown  that 


h  ^2  ^3 


s  Q  s  Q  s>  **^"**2*  3 


n(2/,+i)  “1  “2“  “3 

i«i  *^1,  "'^2  *^3 


.  exp(-l5.  r  ) 

“l^V“3  (12) 

The  coaqiilete  set  of  saiiq>led  values  (occupation  numbers)  is  therefore 

. 

fully  determined  by  the  spectral  components  in  k  _  with 

n^,n2,n^  ranging  in  magnitude  from  zero  to  ^It  follows 

that 


mj^,m2,m  3  '  1  2*'3^nj^,n2,n. 


“r“2,“3 
(13) 


)l  = 
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3*1  THE  DISraXBQIFIQir  OF  TBE  SAMPLED  7ALQEB 


We  now  eonelder  tb«  the»c47namlc  eulMyitnn  coneletlng  of  the  oe- 

cupetlcm  nuMber  9^  ...  /  J  )  ^  ^  “  ^(1  i  i  ^  ‘  ^**'^'*^ 

that  the  velocltiee  of  the  ^  ^  3  ij^nvlduel  ^  ^  ^partlclee  are 


Individual 


statistically  unxvlated  to  these  occupation  nunbers,  ve  conclude  ^lat 
the  probability  of  one  particular  of  the  states  corresponding  to  a 
sequence  of  occupation  nuebers  is  given  by  a  Olbbs  distribution  function 

erp  [-W(..Nj^..,..nj^..)/KP3  (!*♦) 

The  pentutabiUty  of  these  states  is  given  by 


(v)‘- 


tv)'- 


“l'- . “SdiV,)' 


The  probability  density  of  a  sequence  of  occupation  numbers  therefore 
can  be  expressed  as 


P(.n, -  -i— 

n  f .  n  s . 


uy)'-  (n^V)l 

o —  o„.f,  exp(-W/KT)  (l6) 


As  long  as  the  occupation  nxsabers  are  fairly  large  we  can  use  Stirling's 
fomula  for  the  factorials.  [Hie  distribution  then  simplifies  to 


p  eqp(»v/KT}  .  exp 


Frcne  Eq.  (13)  ve  know  bow  the  occupation  nuad>ers  nay  be  eiqpressed  in 
terns  of  Fourier  coxponents.  We  oust  note,  however,  that  Eq.  (13)  has 
twice  as  many  terns  on  the  right-hand  side  as  on  the  left  because  N(1i) 
has  both  real  and  Imaginary  parts.  In  changing  the  varlebles  froa  oc¬ 
cupation  numbers  to  Fourier  conponents  we  therefore  count  only  direc¬ 
tions  of  k  pointing  into  one  henlsidhere  provided  we  want  to  use  n^, 
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tad  M  ladipesdMit  rwiMMB.  Ih  tbl«  mm  and  R^ 
■tdad  tor  tte  iMgInary  and  Um  rMl  pavta  of  M0i),  vlth  alallar  no- 
tntloa  for  tbo  oloetrona.  Baeauao  ttao  Fourlar  caaponanta  aro  llBMKrly 
XiAtttad  to  tha  ooevpotioa  iiurbM*  (aeo  Iq.  12),  tbs  JaeoibUa  of  the  traaa 
fOTMtlon  la  a  eooatant.  Wa  l^refors  obtain  for  the  Joint  distribution 
of  the  real  and  laagalnary  parte  of  the  Paurler  coagoamta 


eacp 


^2  ^3 


-  f  III  ^ 


n^-Og-n^  . 
0  -ilg 


+  (n^  +  n^)  -  2Z(Nj.nj.  +  N^n^)] 


(19) 

Here  we  have  Introduced  (ax^)*^  »  (d|^  ^  «  1 

^  -  6o“/V^ 

D-  la  known  as  the  Debye  length  (Spltaer  1956). 

Ws  laardlately  aee  that  this  la  a  aultldlaienalonal  probability 
density  of  the  Qauasian  type.  We  note  that  the  different  Fourier  conpo- 
nents  enter  through  producta  'of  dlstiributlon  functions  for  each.  There¬ 
fore,  we  conclude  that  the  Fourier  congponents  corresponding  to  different 
wave  nuaibers  are  statistically  independUmt.  In  each  of  the  elenentary 
distribution  functions  we  see  that  even  the  real  and  laaglnary  parts 
corresponding  to  the  saan  wave  nuaber  are  Independent.  We  can  there¬ 
fore  write  down  the  egression  for  the  distribution  of  the  real  parts 
of  R(tc)  and  n(k)  for  one  particular  wave  nuaber  escqparately: 


[R^Z(1  +  2X^Z)  +  n2(l  + 
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\m(t)  -  n(t)l^ 


1 

1 


KT 


23^(1  +  Z) 

1  +2X^(1+  Z) 


(2U) 


This  reduces  to  the  resvtlt  of  Pines  and  Bohm  (19^2)  If  ve  put  Z  «  0, 
l.e.,  if  ure  let  the  Ions  become  smeared  out  to  a  uniform  background 
of  positive  charge. 

We  note  that  the  above  results  are  Independent  of  the  presence 
of  a  magnetic  field  or  tvo-body  collisions  (see  Eq.  10).  This  has 
only  been  shovn  here  for  the  case  of  a  plasma  vith  Coulomb  interaction. 
Similar  conclusions  have  been  reached  by  Dougherty  and  Farley  through 
entirely  different  considerations. 


IV.  THE  SPECTRAL  DISTRIBUTION  OF  THE  SPATIAL  C(»IFONEIITS 

We  next  desire  to  know  how  the  density  fluctiiations  vary  with  time. 

A  relationship  describing  these  time  variations  is  given  by  the  Boltsmann 
transfer  equation.  This  equation  has  been  solved  for  the  case  of  a  plasma 
in  an  external  meignetlc  field  (Chross  1951/  Gordeyev  1952  and  Bernstein 
1958)  •  It  is,  therefore,  here  necessary  only  to  tBdte  over  these  re¬ 
sults  and  apply  them  to  our  problem  directly.  We  do  have  to  supply 
some  sort  of  initial  conditions  to  obtain  the  answer  we  want.  As  we 
shall  see,  these  can  be  found  from  the  results  of  the  previous  section. 

For  completeness  the  method  for  solving  the  Boltanann  equation  is 
briefly  outlined. 

4.1  SOLUTION  OF  TEE  BOLT23(ANN  ^UATION 


The  Boltzmann  equation  for  the  peurtlcle  density  in  phase  space 
is  given  by: 


[l(r,t) 


coll 


-  10 


(25) 


Here  s  >  e/m  for  electrons  and  Ze/M  for  the  Ions.  Assuming  that 
the  deviations  from  a  maxwelllan  velocity  distribution  are  fairly  small, 
we  put 


n(f;V,t)  =  n^(v)  II  +  nj^(r,v,t)]  (26) 

with  n^(v)  a  maxwellian  and  n^^  small  compeured  with  mity.  Again, 
introducing  spatial  Fourier  tmisforms  by 

n,(r,v,t)  =^n  (t,v,t)  exp(-l?r)  (2?) 

t 

and  Laplace  transforms  through 

00 

n^(t,v,s)  n^(^,v,t)  exp(-st)  dt  (28) 

0 


the  linearized  eqxiatlon  becomes 


snj^(iK,v,s)  -n^($,v)  -l5vn^(iJ,v,8)  +  |i 


|_n  (v)  dv  ov 


=  -  vnj^(k,v,8) 


(29) 


The  quantity  n^(^,v)  is  the  initial  value  obtained  by  putting  t  =  0 
in  nj^(^,v,t).  In  the  above,  the  collisions  have  been  tedten  into  ac¬ 
count  by  meeuis  of  a  relaxation  term.  The  quantity  y  is  an  effective 
collision  frequency.  By  combining  the  relaxation  term  with  the  first 
term  on  the  left-hand  side  we  obtain 


sni(lt;V,8)  -►  (s  +  v)  n^(J,v,s)  =  8’n^(^,v,s) 


We  can  therefore  neglect  the  collisions  in  the  calculations  provided 
we  remember  that  wherever  s  appears  e:q)llcltly  it  actually  stands 
for  the  sum  of  s  and  the  collision  frequency. 

To  solve  the  differential  equation  in  v  one  may  Introduce  cylin¬ 
drical  coordinates  with  axis  along  the  magnetic  field  and  with  propagation 
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vector  In  the  1-3  plane,  as  shown  In  Fig.  1.  With  these  w  coordinates 
Eq..  (29)  becomes 


iX(U  COB  ^  +  w  sdn  e  cos  0)]  = 


_  1 
PB 


(30) 


The  Integration  of  Eq.  (30),  which  has  been  performed  by  Bernstein  (1958), 
Is  briefly  outlined  In  Appendix  A. 


3 


FIG.  1.  COORDIKATE  SYSTEM  FOR  SOLVING  EQ.  (29) 


The  solution  Is; 


1  u  dn ' 

n  (h,v,s)=—  \  —  ^  l(ii,s)  -n^(^,v')  . 


60' 


exp|—  [(s  -  Iku  cos  e)i0-0')  -  Dew  sin  e  (sin  0-sln  0 


(31) 


■") 
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If  this  Is  written  out  explicitly  for  electrons  and  ions,  one  obtains 


»  -  ^Oe(0>0')  l2H(^,s)-n(^,s)]  [ 

“o  3 


and 


N^(1c,'^,s)  =  -  ^0^(0, 0*)  +  Z - 2  ^^'[ZN(ic,s)-n(^,B)] 


d0’ 

(32) 

dj?* 


(33) 

where  the  integrating  factors  Ge  and  Gj  follow  from  Bq.  (31). 

We  next  remove  fluctuations  in  velocity  space  by  averaging  over 
all  velocities  by  means  of  the  maxwellian  distribution,  that  is,  we 
form 


n(5,8)  =  /  n^Cv)  n^(^,v,8)  d(v) 

K(t8)  =  /  N^(^)  N^(^;^,8)  d(^)  (3»*) 


This  le8d.8  to  the  eqiuitions 


n(J,s)  =  y^(?,s)  -  ^  {ZN(^,8)  -  n{lt,B)) 


N(]^,s)  =  Yi(^,s)  +  ^  2x5?i  {ZN(?,s)  -  n(^,s)) 


“o 


(35) 


The  complete  expressions  for  Y  and  R  Introduced  both  for  ions  and 
for  electrons  are  given  in  Appendix  B.  Solving  the  set  of  equations 
(35)  one  obtains  .  o  o  i  o 

-  r  23CZ^i)  -  aCZB. 

n(I,s) - 2_i - 2_I_  (36) 


^  -  f 
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(37) 


N(t,e) 


-i.f 

o  0 

1  -  ~  Z^i) 

o  ^ 


From  the  above  two  equatloas  the  variation  of  electron  density  and  Ion 
density  with  time  can  be  foimd  by  taking  the  inverse  Laplace  transform. 
Only  the  time  variation  of  the  electron  density  will  be  computed  in  de¬ 
tail  here,  as  the  procedure  for  dealing  with  the  ion  density  variations 
is  similar. 

For  the  electron  density  the  inverse  transform  becomes 

+i<»fe 

=5^  ^  n(^,8)  e3Q(8t)  ds  (38) 

-!«»+£ 


The  electron  density  at  time  t  can  therefore  be  found,  provided  the 
initled.  conditions  at  time  t  >  0  are  known*  Because  of  the  statis¬ 
tical  nature  of  the  problem,  however,  we  cannot  fix  the  Initleil  condi¬ 
tions.  We  mtist  therefore  resort  to  a  statistical  description.  Hence 
we  form  the  following  expression 

+ii»fe 

<n*(^,0)n(k,t)>^^  =  ~  <n*(]^,0)  n(^,8)>^^e3q)(st)  ds 

-ioo+e 


(39) 

The  spectiTum  is  the  Fourier  transform  of  eui  autocorrelation  function, 
defined  by  (Landau  and  Llfshitz  1958) 

■|  <n*(l^,0)  n(^,t)  +  n(k,0)  n*(l,t)>^^  (40) 

Because  of  the  symmetry  properties  of  n(k,t)  this  turns  out  to  be 
equal  to 


2Re<n*(?,0)n(t,t)> 

bWt 


-  II+  - 


By  means  of  this  e]q>rea8lon  and  by  using  the  Wiener-Khlnchlne  theorem^ 
we  conclude  that  the  spectral  distribution  is  given  by 


—  lln  Re  <n*(^,0)  n(5,8)> 
”  Re(s>^ 


(kl) 


Equation  (4l)  is  the  formal  solution  of  our  problem.  In  order  to  be 
able  to  compute  the  actual  spectra  we  must  study  in  some  detail  the 
function  <n*(J,0)n(^,s)>  \rtiich  is  given  by 


<n*(^,0)y,>^^  (1  .  ^  23^ZR, 

■ _ o  o 

1  -  2X^  (Re  +  Z^i) 

o  ^ 

(h2) 

In  particular,  the  eiqKresslons  for  the  averages  contain  terms  of  the 
form 


and 


<n*(^,0)n^(^,v)>^^ 

<n*({,0)N^(ic,^)>^^  (43) 


In  the  previous  section  we  assumed  that  the  spatial  density  fluctuations 
are  Independent  of  the  velocities  of  the  individual  particles.  As  long 
as  this  is  true,  we  conclude  that  the  ejqpresslons  (43)  are  both  indepen¬ 
dent  of  velocity.  Using  Eqs.  (2l)  and  (22)  we  obtsdn 


<n*(5,0)n^(t,v)  ~  <|n(J)|^>^^ 

0 


,  1+23^Z 

1  _P 

^  1  +  2X^(1+Z) 


<n*(l,0)H^(^,^)>^^  =  <n*(^)N(lt)>^^ 

0 


1 

^  1  +  23^(l+z) 


(44) 
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Substitution  of  these  into  eaqjressions  (lt2)  and  (Ul),  together  with 
the  results  of  Appendix  B  we  finally  obtain  for  the  spectrum 


,  .  n„  Im(-Fe)|l+2X?ZFil2+4X^ZIta(-F£)lFj 

<ln(.,»,  1^.. .  - 


The  function  P,  as  shown  in  Appendix  is  defined  in  the  following 
manner  for  electrons; 


00 

^i  |-  +  expJ-iy(|-)-3i^^^sin^e(l-cos  y) 


1 

+  •jy 


■^COS^©]| 


dy 


Here  we  have  defined  the  quantities 

2 

^  0) 

^  “Scf 


X,  and  Ag  by 


and 


V 

e 


n 


a 


is  the  gyrofrequency  of  the  electrons,  eB/m,  and  Vg  is  the  effec¬ 
tive  collision  frequency  of  the  electrons.  See  Bq.  (29)  and  comments. 
The  expression  for  Fj  may  be  obtained  by  substituting  kX  for  X 
and  for  in  the  expression  for  F^.  k  Is  the  square  root 

of  the  ratio  of  ion  and  electron  masses.  One  must  also  substitute 
the  approp'^iate  collision  frequency  for  ions. 

The  corresponding  expression  for  the  ion  density  fluctuations 
become 


n^  Im(-F^.)|l-^23^Fj^  -f  lm(-F,)4X^z|F 

tfZVo)  ll+23d(F,+ZPi)l^ 

(46) 


and  for  the  charge  density  variations 
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2 

n^e  lin(-Pj)  +  Z  Im(-Fi ) 
nVco  l  l  +  2X^(Fe+ZFi  )1^ 


m 


In  the  next  section  we  go  on  to  discuss  the  properties  of  the 
electron  density  spectra  under  conditions  of  interest  in  ionospheric 
applications . 


V.  DETERMINAriON  OF  ELECTRON  DENSnT  SPECTRA 

Before  starting  to  discuss  the  shape  of  the  spectra  of  electron 
density  fluctuations  for  application  to  scattering  in  the  ionosphere, 
we  must  determine  from  ionospheric  data  the  ranges  of  the  various  param¬ 
eters.  It  turns  out  that  the  spectra  are  quite  straightforward  to  com- 
pute  when  k  is  parallel  to  the  magnetic  field.  This  longitudinal 
ceuse  is  identical  with  the  nonmagnetic  case.  When  ^  is  no  longer 
parallel  to  the  magnetic  field,  it  is  still  possible  to  work  out  exact 
e^qjressions  for  both  F^  and  F.  but  these  involve  slowly  convering 
infinite  sums  with  terms  proportioneLL  to  Bessel  functions  of  imaginary 
arguments  and  are,  in  general,  difficult  to  handle.  However,  in  the 
limiting  cases  of  radius  of  gyration  of  the  charged  particle  small  or 
large  in  comparison  with  k”^,  more  tractable  approximate  expressions 
can  be  used.  In  the  ionosphere,  within  a  fairly  wide  rauiio  frequency 
range,  the  electrons  will  have  a  smaLIl  raidius  of  gyration,  and  the  ions 
a  large  radiuS'  Cf  gyration,  compared  with  k“^. 

For  strictly  transverse  propagation  the  spectrum  cannot  be  worked 
out  without  retaining  some  sort  of  damping.  As  the  extent  of  this  damp¬ 
ing  is  not  known,  the  spectra  can  be  discussed  only  qualitatively  in 
this  case. 

^.1  RANGE  OF  THE  PARAMETERS 

We  shall  assume  that  over  the  height  range  of  Interest  at  present, 
say  from  about  300  km  to  roughly  2000  km,  the  temperature  is  of  the 
order  of  1000°  to  2000°K  (Chapman,  i960).  In  the  actuaJ.  calculations 
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of  the  spectra  we  shall  consider  only  the  case  of  thermal  equilibrium 
between  ions  and  electrons.  As  a  representative  value  for  the  gyro- 
frequency  of  the  electrons  we  shall  take  1  Mc/s.  The  type  of  ions 
present  is  not  at  all  well  known.  At  greater  heights  there  is  some 
reason  to  believe  that  protons  may  predominate^  whereas  at  lower  heights 
oxygen,  nitrogen  and  others  maybe  present.  To  get  some  insight  into 
the  effect  of  the  mass  of  the  ions,  spectra  sire  given  both  for  protons 
and  for  oxygen  ions. 

The  plasma  frequencies  are  taken  to  range  from  0.5  to  5*0  Mc/s 
and  the  radio  frequencies  are  assimied  to  lie  in  the  range  from  50  to 
3000  Mc/s.  One  then  obtains: 


Range  of 

Range  of 


«  "  2KT 


CVi 

10"^  7.10- 

2 

(r)  ’ 

4.10”®  -►  3-10' 

Only  singly  ionized  ions  will  be  considered.  The  values  for 
K  =  (w/m)^  are; 


Hydrogen  :  c  =  43 
Oxygen  :  k  =  172 


5.2  LONGITUDINAL  CASE 

This  case  is  obtained  by  putting  9-^0  in  the  integrals  for 
and  F^  .  First  we  would  like  to  know  whether  collisions  have  to 
be  tedcen  into  account.  By  studying  the  Integrals  we  conclude  that  col¬ 
lisions  are  unimportant  provided; 

=  ^ce  «  1 

and 

X  A .  =  X  .  «  1 

J  Cl 

These  conditions  amount  to  assuming  that  the  mean  free  paths  of  the 
electrons  and  the  ions  are  larger  than  the  scale  of  the  spktial  Fourier 
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component  under  consideration,  i.e,,  larger  than  k""^.  Since  the  mean 
free  path  at  200  km  is  of  the  order  of  several  hundred  metres,  and  be¬ 
cause  the  vavelength  of  the  radio  waves  will  be  at  most  a  few  metres, 
we  conclude  that  collisions  are  unlnqportant.  With  no  collisions  we  find: 

Re(F^)  =  1  -  2X  e:q>(-X^)  ^ 

0 

Im(-Fg)  *=  it^X  e3ip(-X^)  (48) 

Also  we-aote  that"~F^  (X)  =  F^(/cX|. 

In  Eq.  (45)  the  quantity  appears  as  a  factor  and  is  taken 

outside  the  fraction  and  combined  with  the  first  factor.  Dividing 
Eq.  (45)  by  the  resulting  new  first  factor  we  obtain 


kV  ,a>rKT,i  ,1  ,12  ekp(-)t^)ll+24F/,l.X^  «*p(- 

HW  =—  , - ; - 4 - 4-- - 

^  U+2^(F. 


m 


av 


(^^9) 

The  right-hand  side  of  this  was  confuted  on  an  electronic  computer  for 
values  of  X^  ranging  from  300  to  0.03  for  /c  =  43  Md  jc  =  172.  The 
results  axe  shown  in  Fig.  2. 

Let  us  try  to  interpret  these  curves  in  as  simple  physical  terms 
as  possible.  Consider  density  variations  of  one  scale  k"^  only.  IHie 
time  variations  of  the  density  fluctuation  at  this  scale  may  be  thought 
of  as  a  supezposltlon  of  highly  damped,  plane  longitudinal  electron 
waves  with  a  wide  range  of  different  phase  velocities.  These  waves  are 
excited  through  interaction  with  the  microscopic  motion  in  the  plasma. 
For  large  scedes  k’^  the  electrons  are  tied  to  the  ion  motion  and  the 
spectrm  is  identical  to  that  of  the  ion  motion.  This  case  is  repre¬ 
sented  in  Fig.  2  by  the  curve  for  X^  =  300.  As  the  scale  k"^  is  de¬ 
creasing  and  becomes  of  the  order  of  the  Debye  length  the  electron  mo¬ 
tion  is  no  longer  completely  tied  to  the  motion  of  the  ions.  In  ftict, 
the  fast  electron  waves  are  actually  becoming  Independent  of  the  ion 
motion  before  the  slow  ones.  This  explains  why,  in  the  transition  re¬ 
gion,  the  spectral  density  Increases  with  velocity  in  the  velocity  range 
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FIG.  2.  H(X)  -  iil  (12!*I)X 

"o  ■ 

FOR  THE  LONGITUDINAL  CASE. 


<|n(k,a))  PLOTTED  AGAINST  X 
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vhere  the  fast  electron  waves  are  able  to  move  Independently  of  the  ion 
motion.  In  particular,  the  pronounced  dip  in  the  curves  for  s-  1.0 
and  3.0  can  be  explained  this  way.  When  3^  becomes  small  the  spectrum 
tends  toward  a  Oausslan,  as  it  sho\xld  when  the  electrons  become  completely 
free  to  move. 

For  large  Xp  there  is  a  very  sheurp  peak  in  the  curves  for  X  X^. 
This  corresponds  to  the  familiar  electrostatic  waves  at  the  plasma  fre¬ 
quency.  The  contribution  to  the  total  fluctuation  from  this  peedi  is, 
in  general,  negligible. 

It  should  also  be  noted  that  the  widths  of  the  curves  for  leurge 

X  and  for  small  X  are  in  the  ratio  x:”^. 

P  P 

5.3  GENERAL  CASE 


Useful  approximate  expressions  for  F.  and  F,-  may  be  developed 
in  the  general  case  of  d  ^  Let  us  assume  that  collisions  are  un¬ 
important  and  that  we  avoid  the  case  of  strictly  transverse  propagation. 
Then  the  Integral  in  F^  and  .;T  •  can  be  solved  by  eaqptodlng 

V  l> 


sln^  e 


exp 


.2X 


cos  y 


in  a  Fourier  series  and  Integrating  term  by  term. 

The  e^qpresslons  for  the  reeLL  suid  imaginary  parts  of  F.  then  become 

9 


Re(Fj  =  1  - 


2X 


cos  6 


exp 


2X  ,  /itlco 


I  /■«( 


lia(-F^) 


(50) 

where  I^  is  the  Bessel  function  of  imaginary  argument  of  order  n. 
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Because  these  sums  are  rather  difficult  to  handle  we  shall  here 
deal  only  with  approximations  for  leurg^  and  small  radii  of  gyration. 

In  ionospheric  applications,  the  former  case  applies  to  the  ions  and 
the  latter  to  the  electrons,  for  a  fairly  large  range  of  wave  vectors 
When  the  radius  of  gyration  is  small,  >  1  and  only  the  zeroth 
order  terms  in  the  sums  above  need  be  considered  as  long  as  X  ^  also 
exceeds  X.  Hence: 

X 

cos  6 

^  “I'  >  5 

0 

(51) 

Note  that  these  equations  could  have  been  obtained  from  the  longitudinal 
case  by  replacing  X  by  x/cos  9. 

When  the  radius  of  gyration  is  large  an  approximation  can  best 
be  obtained  from  the  original  expression  for  rather  than  from 
Eqs.  iy>).  In  this  case  X^  <  1,  and  the  expression 

[sln^e  ”1 

- 5-  (1  '  cos  y) 

2x;  J 

will  be  appreciably  different  from  zero  only  when  y  =  2jtn,  provided 
sin  e  >  X^.  The  contributions  around  y  =  2rtn  are  obtained  by  e]q)and- 
Ing  1  -  cos  y  to  second  order,  i.e.,  by  putting 

1  -  COB  y  «  i  (y  -  2jra)^ 

The  expressions  for  the  real  and  Imaglneo'y  parts  of  F,  specialized 
for  the  ions  through  the  substitution  (see  p.  l6) 

X  ^  xX 

X  -►  a:"^ 

e  e 

become 
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Re(Fi)  =  1  -  2  /cX  exp  [-  x^X^l 


exp  Bin^e  coB^e] 

A« 


Bin  [2imc^  Bin^e 


Im(-P.)  =  jt^  /cX  exp  [-  /c^X^]  exp  [-(~)^  ein^e  coB^e] 

(52) 


m=l 


2  2  X 

COB  (2jtnvc  sin  0 


The  effect  of  the  magnetic  field  is  to  modulate  the  functions  obtained 
in  the  longitudinal  case.  It  is  only  when  6  is  fairly  close  to  90° 
that  that  this  modulation  becomes  appreciable. 

Because  of  the.  number  of  pBuraneters  Involved  it  is  not  possible 
to  (Knpute  sets  of  curves  '^Id  for  as  many  different  cases  as  for  the 
nonmagnetic  case  of  Fig.  2.  In  order  to  discuss  the  effect  of  the 
magnetic  field  on  the  spectral  distribution^  we  therefore  choose  one 
particular  set  of  parameters  and  study  the  change  of  the  spectrum 
with  angle  Q . 

Vfe  choose  parameters:  ^ 


X^  =  10 
Xp  =30 

K  =1*3  (hydrogen) 


Let  us  discuss  the  significance  of  this  choice. 


As  mentioned  above, 
3 


the  electron  gyro- frequency  is  taken  as  1  Nc/s.  Because 
this  means  that  the  plasma  frequency  is  3  Mc/s,  which  is  at  least 
of  the  rlglit  order  of  magnitude  above  the  peak  of  the  F- layer  of  the 
ionosphere.  If  the  teDq)erature  is  1000°K,  then  (^pp)^  ®  5*7  •  lO’^. 


With  X,  =  10  in  a  backscatter  experiment  the  radio  wavelength  would 
correspond  to  3.5  meters,  which  is  typical  of  the  wavelengths  used  at 
present . 
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Substituting  Bqs.  (5I)  and  (52)  into  Eq.  (45)  and  normalizing  as 
in  the  longitudinal  case  (Eq.  49)  we  can  again  compute 

H(X)  =  ^  (^)^ 
n^  m 

eigalnst  X  on  an  electronic  coo^uter.  The  computation  was  carried 
through  for  6  =  8o°,  85°  and  88°,  and  the  "ionic"  parte  of  the  spectra 
are  shown  in  Fig.  3*  The  longltudinsd  case  for  =  30  is  shown 
dotted. 

We  see  that  as  6  approaches  90°  the  spectrum  develops  toward  a 
line  spectrum  with  peaks  at  the  gyro-frequencies  of  the  ions,  and  peaks 
at  all  the  harmonics  of  this  frequency.  Because  of  the  computer  time 
required  the  range  of  the  power  density  was  limited.  The  frequency 
range  was  also  reduced,  for  the  same  reason,  to  the  reuge  of  X  where 
the  spectrum  is  essentially  of  ionic  nature. 

The  spectrum  for  6  =  85°  is  displayed  on  a  linear-linear  scale 
in  Fig.  4. 


FIG.  4.  PLOT  OF  H(X)  AGAINST  X  ON  A  LINEAR  SCALE. 
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VI.  DISCUSSION  AND  CONaUSION 


The  present  work  has  shown  how  the  spatial  Foiurler  conponents  of 
electron  density  fliictuatlone  and  their  spectral  distribution  ms^  be 
cooputed  for  a  plasaa  in  thermal  equilibrium  in  an  external  magnetic 
field  when  Coulomb  Intex^ctlon  is  taken  into  account. 

In  the  ionosphere  we  will  be  concerned  with  only  Z  s  1.  Thus, 
if  the  "scale"  is  leurger  than  the  Debye  length,  the  total  fluctua¬ 
tion  is  found  to  be  half  that  found  in  the  absence  of  interactions. 

If  the  "scale"  is  smaller  than  the  Debye  length,  the  toted,  fluctuation 
is  not  Influenced  by  Interactions.  The  total  fluctuation,  and  hence 
the  total  scattered  power,  is  independent  of  the  presence  of  a  stav.ic 
magnetic  field. 

For  physical  conditions  which  are  believed  to  prsvall  in  the  iono¬ 
sphere,  and  for  wavelengths  (scales)  of  the  order  of  a  few  meters,  the 
spectral  distribution  is  not  influenced  by  the  presence  of  a  magnetic 
field  unless  the  wave  vector  ^  is  nearly  perpendicular  to  the  mag¬ 
netic  field.  As  ^  approaches  perpendicularity,  sharp  resonance 
peaks  are  developing  rapidly  round  the  ionic  gyro- frequency  and  multiples 
thereof. 

It  therefore  appears  that  it  will  be  possible  to  study  the  ionic 
constituents  in  the  ionosphere  by  carefully  designed  radio-wave- scatter 
esqperlments .  The  temperature  can  also  be  found  by  stxidylng  the  spectral 
width,  either  of  the  "ionic"  spectrum  or  of  the  "electronic"  spectrum. 

The  electron  density  may  be  found  from  the  total  scattered  power. 

The  possibility  of  ewtually  detecting  the  peaks  at  the  ionic 
gyro- frequencies  depends  on  the  frequency  resolution  available.  In 
normal  backscatter  obseirvations  with  relatively  short  pulses  there  is  a 
chance  that  the  details  may  be  washed  out  because  of  the  width  of  the 
spectrum  of .the  transmitted  wave.  It  therefore  appears  that  some  sort 
of  blstatlc  continuous-wave  e]q;>erlment  ought  to  be  looked  into  and 
'thoroughly  discussed. 

When  several  ionic  components  are  present  the  above  theory  can  be 
extended  without  any  difficulties  in  principle,  but  the  work  becomes  more 
laborious  and  is  omitted  here. 
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APFERDIX  A 


13le  solution  of  the  hcmogeneous  part  of  Eii.  (30)  is: 


Cq(W/U) 


exp 


1 

—  [(s  -  Iku  cos  6)0  -  ikw  sin  6  sin  0] 
HB 


i 


(Al) 


Instead;  of  C^(v,u)  we  substitute  a  function  C(w,u,^}  and  try  to  solve 
the  Inhomogeneous  eq.uatlon.  It  Is  then  found  that: 


C(w,u,0) 


^flied 

limit 


1 

—  [(s  -  Iku  cos  6)0* 


-  Ikw  sin  6  sin 


-  n-^(k,v' ) 


I  d0' 


(A2) 


nie  full  solution  therefore  becoMs: 


n,(w,u,0)  -  —  \  expi—  [(s  -  iku  coo  e){0  -  0*)  -  .'a 

"  tiL  r 


tv  sin  6 


limit 


(sin  0  -  sin  0*)] 


^a* 

^B(k,s) 

.dv' 


n  (k,v')(>  ^0 


(A3) 


where  v  is  obtained  from  v  by  putting  0^0'.  Because  physical  con^ 
slderatlons  require  that  the  function  be  single -valued,  we  must  construct 
a  solution  which  is  periodic  In  0.  The  only  possible  solution  Is  then 
the  one  presented  In  Eq.  (31)* 
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APFKRDIX  B 


Prom  Eqs.  (32),  (33)  and  (35)  we  see  that; 

Y,  -  -  y  ^  )nQ(v)nj^(k,v' )d( v)d0» 


-If  k7'Ge(0,0' 


)n^(v)d(v)d0' 


(Bl) 


and: 


Y.  -  . 


«i 


y  ^  G.(0,0*)N^(V)N^(k,V' 


)d(V)d0' 


Y  ao 


-if  kV'G.(0,0' 


)NQ(V)d(V)d0' 


(B2) 


The  expressions  for  R  may  he  Integrated  imawdlately,  but  not  so  those 
for  Y  because  of  the  stochastic  nature  of  n  (k,v')  emd  ^(k^V').  Before 
integrating  we  must  form  the  required  averages  <n*(k)Y,>  and 

BV 


<n*(k)Yi>^^  and  make  \ise  of  Eqs.  (44).  The  following  results  are  then 
obtained: 


<n*(k)Ye> 


1  l-t-2X^Z 

_  _ E_ 


av 


V  1  +  23^(1  +  Z) 
P 


y  f  G,(0,0' 


)n^(v)d(v)d0' 


(B3) 


<n*(k)Yi> 


1  2X^Z 

JL 


V  1  +  2X^(1  +  Z) 
P 


I  1 (<»-«»' 


)N^(V)d(V)d0' 


(B4) 


Bernstein  (19^)  has  shown  how  this  type  of  Integral  can  be  evaluated. 
Because  the  procedure  is  fairly  straight  forward  the  details  are  omitted 
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here. 


!nie  IntegreLL  In  (B3)  becomes: 


o  2  2 

[sin  0(1  -  cos  3)  +  cos  0] 


EBt®  ^ 


The  Integral  In  (B4)  Is  of  exactly  the  sane  form,  ve  have  only  to  sub¬ 
stitute  parameters  applying  to  Ions  instead  of  electrons.  The  Integral, 
which  we  will  tenn  a  Gordeyev  Integral,  will  now  be  denoted  by  g(d,s/n), 
and  we  put  indices  <  or  i  to  Indicate  whether  it  applies  to  electrons 
or  ions. 

In  terms  of  this  Integral  the  expressions  (B3)  and  (Bh)  became: 


<n*(k)Yi>^v 


n  1  +  23rz 
_2 _ P 

fleV  1  +  2X^(1  +  Z) 
n  1  +  23^2 

_2 _ E _ 

n^V  1  +  2X^(1  +  Z) 


(B5) 

(B6) 


The  Integrals  occurring  in  R  ,  and  R  ^  nay  also  be  expressed  in  terms  of 
the  Gordeyev  integral: 


Our  final  formulae  (4^),  (46)  and  (4?)  follow  from  these  by  simple 
substitution. 


